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Abstract 

The Ks — ► 77 decay rate has been measured with the NA48 detector using a 
high intensity short neutral beam from the CERN SPS. The measured branching 
ratio BR{K S -► 77) = (2.78 ± 0.06 stat ± 0.04 SJ/st ) x 10~ 6 , obtained from 7461±172 
Ks 77 events, is significantly higher than the 0(p 4 ) prediction of Chiral Per- 
turbation Theory Using a Kl beam the ratio p^^tt^I ) = (2.81 ± 0.01 stat ± 
0.02 SJ/st ) x 10 -3 has been measured. 
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1 Introduction 

The decays Ks,l ~* 77 an d Ks,l ~~ > 7r°77 are important probes of Chiral Pertur- 
bation Theory (xPT), an effective field theory of the Standard Model at low energies. 
The decay amplitude of K$ — > 77 can be calculated unambiguously at the leading 0(p 4 ) 
order of the perturbative expansion giving a branching ratio of 2.1 x with an un- 
certainty of only a few per cent |1]|q The theoretical predictions for Kl — > 77 are much 
less accurate |§. The K$ — > 77 decay rate has hitherto been considered as an exam- 
ple of good agreement between the xPT prediction and the experimental measurement 
(2.4±0.9) x 1(T 6 , obtained by the NA31 experiment Q. This result was recently improved 
by NA48 to (2.58 ± 0.42) x 1CT 6 with data from a two-day test run with high inten- 
sity short neutral beam, collected in September 1999. However, the previous experimental 
results are not precise enough to resolve the higher order effects which are predicted by 
the xPT to be at most of the order of m| r /(47rF 7r ) 2 ~ 20% of the 0(p 4 ) decay ampli- 
tude [|Xj| . Using data taken in a 40 day run with the high intensity short neutral beam in 
the year 2000, for the first time an accuracy of a few % has been achieved. A comparison 
of the experimental result with the 0(p 4 ) X-PT prediction allows the estimation of the 
0(p 6 ) loop contributions. According to X-PT, these could also contribute to 0{jP) effects 
in the Kl — > 7r°77 decay which up to now are believed to be dominated by vector meson 
exchange. 

2 Principle of the measurement 

The NA48 experiment has been designed to measure direct CP violation in neutral 
kaon decays and comprises Kg and Kl beams created at two targets, near and far from 
the decay region, which can be run together or separately || . 

In previous measurements of Ks — > 77, the decay region used contained almost 
twice as many decays coming from Kl — > 77. Their subtraction using the Kl flux es- 
timated from e.g. Kl — > 7r°7r 7r° decays introduces an uncertainty of more than 4% due 
to the poorly known Kl — > 77 branching ratio. In the measurement presented here, this 
problem is overcome by using a novel method based on the observed ratio of Kl — > 77 to 
Kl — > 7r 7r°7r°, as a function of energy and decay position, in the Kl beam (far-target) run 
performed immediately prior to the Ks beam (near-target) data taking. As the near- and 
far-target beams are almost collinear and have similar momentum spectra, the measure- 
ment of Kl — > 7T°7rV decays in the near-target beam, as a function of energy and decay 
position, and the ratio A^" r ^ 77 /A^^ 0ff07r o measured from the far-target beam allow the 
direct evaluation of the number of Kl — > 77 decays in the near-target beam. By using 
the same ranges of energy and decay position only a very small acceptance correction is 
needed due to the slight differences in geometries of the two beams. A by-product of this 



method is a precise measurement of 



r(jr L -> 77 ) 



3 Experimental set-up and data taking 

The near-target beam contains K° mesons and hyperons. It was produced on a 
2 mm diameter, 400 mm long beryllium target by 400 GeV protons at a production 
angle of 3.0 mrad in the vertical plane. In this special high intensity near-target mode, 
the primary proton beam was attenuated and collimated to the desired intensity, far 
upstream of the target. The target is followed by a sweeping magnet packed with tungsten- 
alloy inserts in which protons not interacting in the target are absorbed, and by a 0.36 

^ These results are similar to those obtained in earlier calculations based on a phenomenological ap- 
proach H 
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cm diameter collimator. In the usual NA48 experimental configuration the 1.5 m long 
collimator is followed successively by an anti counter (AKS), by an 89 m long evacuated 
tube terminated by a thin Kevlar window and by a helium filled tank which contains the 
drift chambers and a spectrometer magnet. However for these data the AKS counter, the 
Kevlar window and the drift chambers were removed and the helium tank was evacuated. 

The present analysis is based on data recorded during a 40-day run in 2000 with 
a beam intensity of about ~ 10 10 protons hitting the target during the 3.2 s long SPS 
spill. This is a factor of ~ 300 higher than the usual near-target beam. This near-target 
run was preceded by a 30-day run with the far-target beam alone and the same detector 
configuration^], devoted mainly to studies of systematic effects in the direct CP- violation 
measurement |7], || . 

The detector elements used in the analysis are the following: 

— A liquid krypton calorimeter (LKr) || is used to measure the energy, position and 
time of electro- magnetic showers initiated by photons (7). Around 20 t of liquid 
krypton at 121 K are used as an ionization detector. The calorimeter has a structure 
of 13212 square towers of 2 x 2 cm 2 cross-section and 127 cm length (27 radiation 
lengths) each. The cells are formed by copper-beryllium ribbons, 1.8 cm wide and 
40 fim thick, stretched longitudinally. The ionization signal from each of the cells is 
amplified, shaped, and digitised by 10-bit FADCs at 40 MHz sampling frequency. 
The dynamic range of the digitisers is increased by gain- switching amplifiers which 
change the amplification factor depending on the pulse height. The energy resolution 
is 

a(E)/E ~ 0.090/E © 0.032/v 7 ^ © 0.0042, (1) 

where E is in GeV. The read-out system was calibrated by a charge pulse every 
burst during data taking. The final calibration of the energy scale is fixed by fitting 
the effective edge of the collimator seen in the data to be the same as in the Monte 
Carlo. The relative calibrations of the individual cells were determined by using 
K e 3 decays during the 1998 run and checked to be similar in the 2000 run with 
7T°, 77 — > 77 decays produced in thin targets in a special 7r~-beam run. The position 
and time resolutions for a single photon with energy larger than 20 GeV are better 
than 1.3 mm and 300 ps, respectively. 

— A sampling hadron calorimeter composed of 48 steel plates, each 24 mm thick, inter- 
leaved with scintillator planes is used to measure hadronic showers with a readout 
in horizontal and vertical projections. 

— A scintillator hodoscope upstream of the two calorimeters helps to detect charged 
particles. It is composed of two planes segmented in horizontal and vertical strips 
arranged in four quadrants. 

A more complete description of the apparatus can be found elsewhere [BJ . 

The trigger decision, common to 77, tt 0/ k° and 7r 7r°7r° decays, is based on quantities 
which are derived from the projections of the energy deposited in the electro-magnetic 
liquid krypton calorimeter |l(J. The trigger required that the total deposited energy E tot 



is larger than 50 GeV, the radius of the centre of energy is smaller than 15 cm and the 
proper life time of the kaon is less than 9 Kg lifetimes downstream of the collimator. The 
trigger efficiency, (99.8 ±0.1)%, is measured using a minimum-bias sample triggered by a 



2 ) with exception of the proton beam energy which was 450 GeV, production angle which was 4.2 mrad 
and of the spectrometer magnet which, unlike in the near-target run, was powered in order to create 
occupancy conditions as close as possible to those of the direct CP-violation measurement. 
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scintillating fibre hodoscope placed in the liquid krypton calorimeter and is equal for 77, 
7r°7r° and 7r°7r 7r° decays in both the near- and far-target runs. 

4 Event selection 

The energies and positions of electro-magnetic showers initiated by photons, mea- 
sured in the liquid krypton calorimeter, are used to calculate the kaon energy and decay 
vertex. To select K — > 77 candidates, all events with > 2 energy clusters are considered. 
Pairs of clusters which are in time within 5 ns and with no other cluster with energy 
> 1.5 GeV closer in time than 3 ns with respect to the event time are selected. The event 
time is computed from the times of the two most energetic cells of the selected clusters. 
In addition, the 77 pair must pass the following cuts: 

— The energy of each cluster must be greater than 3 GeV and less than 100 GeV. 

— The transverse distance between two clusters is required to be greater than 10 cm. 

— The total energy of the selected cluster pair is required to be less than 170 GeV and 
to be greater than 70 GeV. 

— The radius of the centre of energy, 



Rc = T7e, (2) 

is required to be less than 7 cm, where E iy x iy yi are the i — th cluster energy, x and 

y coordinates at LKr, respectively. 
— The energy deposited in the hadron calorimeter must not exceed 3 GeV in a time 

window of ±15 ns around the event time. 
In order to determine the K$ and Kl fluxes in the beam, the decays K$ — > 7r°7r° and 
K L — > 7r°7T 7r have been used. Similar conditions as for 77 cluster pairs are applied to 
groups of 4 or 6 in-time clusters for 2tc° and 3n° respectively. 

The decay vertex position z vertex of a kaon decaying into photons or 7r°'s is calculated 
using the kaon mass constraint 



V / Ei,j,i>j EjEjKxj ~ Xj) 2 + (yj - yj) 2 } 

^vertex ^Lfer 1<5J 

m K 

where ZLkr is the z coordinate of the LKr calorimeter with respect to the end of the 
collimator. 

In 27t° and 3ir° samples a x 2 cu t of 27 (~ 5<r) is applied to the invariant masses of 
photon pairs calculated using z vertex . For 2tt° decay, the x 2 variable is defined as 



x 2 = 



(mi + m 2 )/2 - m^o 



1 2 



+ 



(mi - m 2 )/2 



1 2 



(4) 



where a± are the resolutions on (m x ± m 2 )/2 measured from the data and parametrised 
as a function of the lowest photon energy. For 7r 7r 7r° decays an equivalent variable is 
defined in the space of three invariant masses. 

Due to the large background from K s — > 7r°7r° decays the study of K s 77 can 
only be carried out in a restricted decay region. The maximum 77 mass from a K s — > 7r°7r° 
decay where two photons escape detection is 458 MeV/c 2 which would lead to an apparent 
vertex shift of 9 m downstream. However, in case some showers overlap and only one or 
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no photon misses the detector this shift can be significantly smaller. Therefore the decay 
region is defined to be between -1 m and 5 m with respect to the exit of the collimator. 

In order to remove events with hadronic or overlapping showers from the 77 sample, 
the shower width is required to be less than 3a above the average value for photon showers 
of a given energy. This cut, which is calibrated from showers in K$ — > decays, removes 
< 0.5% of good K — > 77 events. 

In the near-target data, 19819 events passed the 77 selection cuts while in the far- 
target data more than 6 x 10 4 events were used. Normalisation channels were down-scaled, 
retaining several hundred thousand events for analysis. 

5 Background subtraction 

The 7T°7r and 7r°7r 7r° samples are essentially background free. In the 77 samples 
the following backgrounds have been identified and subtracted: 

— A small hadronic background visible as a tail in the shower width distribution. In 
the near-target sample this background originates from interactions of the beam 
particles in the collimator. In the far-target beam, which is cleaned by three col- 
limators, the background from this source is negligible. However, pairs of charged 
hadrons generated in the far-target beam upstream of or inside the final collimator 
and opened vertically by the last sweeping magnet have been found to contribute 
to the Kl — > 77 sample. 

— Accidental 77 pairs visible as a flat continuum in the distribution of the shower 
time differences. 

— Residual K$ — > 7r°7r° decays with one shower overlap and one low-energy photon 
missing the calorimeter. 

— Dalitz decay K° — > ee'j with an overlapping ee pair or one shower with energy 
below the detection limit of 1.5 GeV. In the far-target run, due to the presence of 
magnetic field, ee pairs do not overlap. 

The last three background sources contribute significantly only to the near-target 77 
sample. 

The hadronic and accidental backgrounds are subtracted using the distribution in 
radius of the centre of energy {Rc), which indicates the deviation of the direction of the 
apparent decay source from the beam axis. A control region is defined with Rc between 
8 cm and 10 cm. In order to obtain the background-extrapolation factor to the signal 
region, Rc < 7 cm, a set of model distributions has been extracted from the data by 
modifying certain selection requirements. 

In the near-target data, in order to study the hadronic background component two 
different samples have been selected. In one, a hit in the hadron calorimeter is required 
and the shower width cut is released (extrapolation factor 1.5). In the other sample, 
as collimator scattering can produce more than two particles, an additional shower in 
the calorimeter is required and the decay region is shortened by 1 m in order to avoid 
contamination from K$ — > 7r°7r° decays, while all other cuts are unchanged (extrapolation 
factor 3.1) (Fig. [I] left). The model of accidental background is obtained by requiring 
the two showers to be out of time by more than 5 ns and less than 10 ns (extrapolation 
factor 2.7). The median value of the extrapolation factors from three different models 
gives an inclusive hadronic and accidental background estimate of (2.1 ± 0.7)% of the 
Ks —>■ 77 sample where the uncertainty corresponds to the spread of the extrapolation 
factors (Fig. [I] right). 
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In the far-target data the background model sample is obtained by requiring a hit 
in the hadron calorimeter and a hit in the scintillator hodoscope giving an extrapolation 
factor of 2.3. Using this factor and a conservative uncertainty leads to a background 
estimate of (0.6 ± 0.3)% (Fig. |2|). In both background estimates the amount of beam halo 
in the control region is taken into account and is estimated from 7r°7r° and samples. 




Figure 1: Left: Example of background subtraction in the near-target 77 sample using as 
a model 77 candidates with an extra shower. Right: Comparison of the centre-of-energy 
distributions of background models in the near-target data. 
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Figure 2: Background subtraction in the far-target 77 sample. 



In order to subtract the Kg — > 7r°7r° background a sample of 1.2 x 10 9 decays 
was generated by Monte Carlo. The Monte Carlo sample was normalised using fully 
reconstructed Kg — > 7r°7r° events and checked by comparing it to the data in the 10 m 
following the decay region where Kg — > tc°tt background largely dominates the 77 sample 
(Fig. [|). The estimated background amounts to (2.1 ± 0.4)% of the Kg — > 77 sample and 
populates mainly the last metre of the decay region. The uncertainty reflects a discrepancy 
between data and Monte Carlo in the simulation of overlapping showers. 

The contribution of K° —>■ eej decays to the 77 sample was estimated using Monte 
Carlo to be 0.4%. In the same way, the numbers of accepted 7r°7r° (0.9%) and 7r 7r 7r° 
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(1.1%) decays with a subsequent 7r° — > eej decay with respect to pure photon final states 
were estimated. The overall correction to the BR(Ks — > 77) due to accepted Dalitz decays 
amounts to (1.5 ± 0.3)%, where the uncertainty is due to the simulation of overlapping 
showers. 




( cm ) Z m (cm) 



Figure 3: The z veriex distribution of the near-target 77 sample with logarithmic (left) 
and linear (right) scales compared to Monte Carlo samples of Ks — > 77, Kl —>■ 77 
and Ks — > 7r°7r° background, normalised to the flux and using the measured values of 
r(^L^°V>i ) anc ^ BR(-Ks — > 77) . The shapes of the distributions before summing are also 
indicated. 



6 Results 

Binned maximum likelihood fits have been used to measure the ratio 
and the K$ — > 77 branching ratio in 10 kaon energy and 6 longitudinal vertex position 
bins. Acceptances were calculated using Monte Carlo with detector simulation based on 
GEANT [[Ll[] where the LKr calorimeter response was simulated using a pregenerated 
shower library. 

The ratio of Kl — > 77 and Kl — > 7r°7r 7r° decay rates, calculated from far-target 
samples and corrected for acceptance and for hadronic background as described in the 
previous section, is 

1 [Kl — > 7T u 7T u 7r u J 

where the systematic uncertainty accounts for the uncertainties in the background sub- 
traction and simulation of the calorimeter response in the acceptance calculation. The 
latter were checked by varying selection cuts, especially the shower-width cut in the 77 
sample and the x 2 cu t i n the 7r 7r°7r° sample, which are specific to each of these decay 
channels. In addition it was checked that the result is insensitive to various parameters, 
e.g. longitudinal vertex position and kaon energy. 

The number of Kl — > 77 events in the near-target 77 sample is given by 

]\rnear p nrnear / R \ 
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Hadronic background and accidentals 
Background from K$ — > tt°tt 
Had. background in far-target beam 
Dalitz decay correction 
Acceptance uncertainty systematic 
Acceptance uncertainty statistical 
Trigger efficiency 


-2.1 ± 0.7 % 
-2.1 ± 0.4 % 
0.9 ± 0.4 % 
1.5 ± 0.3 % 
± 0.3 % 
± 0.6 % 
± 0.1 % 


Total 


-1.8 ± 1.2 % 


Statistical uncertainty data 
BR(K S -> 7r°7T ) 


± 2.0 % 
± 0.9 % 



Table 1: Summary of all corrections and uncertainties in the BK(Ks — > 77) measurement. 



where iV are the background corrected numbers of events from the near and far- 
target runs. Ra is a double ratio of acceptances 

A near A f ar 

_ -n-K L ->'ry 1 /1 AT L ->77 ,j\ 
A A near ' Afar \ I 

which is close to unity (~1.05) because the same decay region was used for the far and 
near-target samples. For this reason only the background subtraction uncertainty in (^) 
enters the calculation of BR,(Ks —>■ 77), while the acceptance calculation uncertainty 
cancels between near and far-target samples. 

The summary of all corrections and uncertainties in the BK(Ks — > 77) calculation 
is given in Table |l|. The acceptance uncertainty reflects the level of knowledge of non- 
Gaussian tails in the energy reconstruction accounted for in the acceptance calculation. 

From 19819 77 candidates, 7461 ± 172 have been estimated to be Ks — > 77 decays, 
giving a ratio of decay rates 

Using BR(K~s — > 7t tt ) taken from |0] one obtains 

BR(K S -► 77) = (2.78 ± 0.06 stat ± 0.03 syst ± 0.02 ex . t ) x lO" 6 (9) 
= (2.78 ± 0.07) x 10" 6 (10) 

where the external uncertainty is that quoted in Table [I] for the Ks — > 7r°7r° branching ra- 
tio. Systematic checks have confirmed that the result is stable with respect to longitudinal 
vertex position, kaon energy, run period and other parameters. 

This result is in agreement with previous measurements, but has much better pre- 
cision. It is the first measurement of BR{Kg 77) to show significant difference with 
respect to the 0(p A ) calculations of xPT and indicates that higher order corrections 
increase the decay rate by about 30%. 
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